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Abstract 

This is a comment on the papers N. D. Elkies, Amer. Math. 
Monthly 110 (2003), 561-573 and Cvijovic and J. Klinowski, J. Corn- 
put. Appl. Math. 142 (2002), 435-439. We provide an explicit ex- 
pression for the kernel of the integral operator introduced in the first 
paper. This explicit expression considerably simplifies the calculation 
of S(n) and enables a simple derivation of Cvijovic and Klinowski's 
integral representation for £(2n + 1). 



1 Introduction 

In the remarkable paper [I] Beukers, Kolk and Calabi showed that for all 
n > 2, 

*<») - £ (3FTTF = ©" *" (L1) 

k=— oo 

where 5 n = Vol(A n ) is the volume of the n-dimensional convex polytope 

A„ = . . . , u n ) : u,i > 0, Ui + u i+ i < 1} . (1.2) 

Here Ui are indexed cyclically modulo n, that is u n+1 = u±. 

In a subsequent paper [3], Elkies gave an elegant method for calculating 
the volume S n (earlier calculations of this type can be found in [8]). If we 
introduce the characteristic function Ki(u,v) of the isosceles right triangle 



1 



{(u, v ) : u, v > 0, w + f < 1} that is 1 inside the triangle and outside of it, 
then [1] 

^ = / ••• / IT^iKj du i ■ ■ - du n = / du ± / du 2 Ki(u 1 ,u 2 ) 
Jo Jo i=1 Jo Jo 

... I du n -i Ki(u n -2, u n -i) / du n Ki( ) i^l(-U n ,Mi). (1.3) 

Jo Jo 

n 

Let us note that Y[ Ki(v>i, u i+ i) is the characteristic function of the polytope 

i=l 

A n and just this property enables us to extend the integration domain from 
A n to the n-dimensional unit hypercube and obtain (II. 3p . At that the value 

n 

of Yl Ki(ui,u i+ i) on the boundary of A n is, in fact, irrelevant and in the 
above integral we can assume 

K x {u,v) = 9{l-u-v), (1.4) 

where 8(x) is the Heaviside step function with #(0) = 1/2 (it doesn't matter, 
as far as the integral (11. 3p is concerned, what value is used for #(0), the choice 
#(0) = 1/2 will be convenient for us in the following). 

Ki(u,v) can be interpreted [1] as the kernel of the linear operator T on 
the Hilbert space L 2 (0, 1), defined as follows 

(£/)(«)=/ K 1 (u,v)f(v)dv= [ f{v)dv. (1.5) 
Jo Jo 

Then (11. 3p shows that 5 n equals just the trace of the operator T n : 

5 n = K^u^m) dm. (1.6) 
Jo 

The kernel K n (u,v) of T n obeys the recurrence relation 

K n (u,v) = K x (u,u 1 )K n ^ 1 {u 1 ,v)du 1 . (1.7) 
Jo 
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2 An explicit expression for the kernel 



No solution of the recurrence relation (11. 7p was given in [J]. It is the aim of 
this short note to provide the solution. When Ki(u, v) = 9(1 — u — v), it is 
given by 



)2n-2 



K 2n (u,v) = 

-^2n-l 



X 



+ 



U + V 



R 



2n-l 



2 

It + V 



and 



^2n+l(M,w) 



X 



+ 



R2n 
R2n 



■1) H 



U + V 



(2n - 1)! 

f -E-2n-l 
+ -^2n-l 



)2n-l 



U — V 

2 

v — u 



6(u-v) 
(v — u) 



{2.1) 



(2n)\ 



2 

u + v 



+ ^2n 



U — V 



E- 



2n 



U + V 



9(1 -u-v 
9(u + v-l 



(2.2) 



In these formulas, R n (x) are the Euler polynomials (see e.g. [7], [2]) and 6*(x) 
is the Heaviside step function with 9(0) = 1/2. 

After they have been guessed, it is quite straightforward to prove by 
induction that (12. ip and (12.21) obey the recurrence relation (11.71) (see the 
Appendix). 

No systematic method was used to get (12.11) and (12.21) . They have indeed 
been guessed. We simply calculated a number of explicit expressions for 
K n (u,v) using the recurrence relation ( 11. 7ft and the expression ( 11.41) . and 
then tried to locate regularities in these expressions. The appearance of 
Euler polynomials suggests that K n (u,v) should have a nice two- variable 
Fourier decomposition and a likely systematic approach could be based on 
this fact [5]. 
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3 Applications of the kernel 

Having (12. ip and (12. 2p at our disposal, it is easy to calculate the integral in 
(11. 6p . Namely, because 

o2n-2 

K 2n (u, u) = (-1)" — [E^u) + £ 2n _i(0)] , (3.1) 

(2n — 1)1 

and 



we get 

$2n 



K 2n (li, U) du = (-1)" (2^,^, ^-1(0) (3.3) 

(note that -E 2 n(0) = -E 2n (l) = 0). But E 2n -i(0) can be expressed through 
the Bernoulli numbers 

E 2n „ l (0) = -^(2 2n -l)B 2m (3.4) 

and using the relation S(2n) = (1 — 2~ 2n ) C(2^), we reproduce the celebrated 
formula for the Riemann zeta function at even-integer arguments 

o2n-l 

a2n) = (-ir +1 —^ n B 2n . (3.5) 
In the case of S(2n + 1), we can use the identity 

2 2n 

K 2n+1 (u, u) + K 2n+1 (l -u,l-u) = (-1)" ^yy ^(1/2), (3.6) 

which follows from ( 12. 2p . and get 

1 

^ J [K 2n+1 (u, u) + K 2n+1 (l -u,l- u)} du 

E 2n (l/2). (3.7) 



>2n+l 





o2n-l 



v ' (2n)! 

Expressing E 2n (l/2) in terms of the Euler numbers 



£ 2n (l/2) = 2- 2n £ : 



2/i • 
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we get from (I3.7p Euler's other celebrated formula (note that formulas (9) 
and (10) in |3J, analogs of our ( 13 .5p and ( 13. 8p . contain typos) 



5(2n + l) = (-1)" 



2(2n) 



7T 



2n+l 



2n- 



(3.8) 



Interestingly, our expression for K n (u,v) allows us to re-derive Cvijovic and 
Klinowski's integral representation [3] for ((2n + 1). We begin with the 
formula (for more details, see [TTj ) 



C(2n+1) 



22n+l | 



2 2n+l _ l 2n 



In {x x ■ ■ ■ x 2n 



u 2n 



1 — x{ ■ ■ ■ X 



(3.9) 



2n 



where n 2n is the 2n-dimensional unit hypercube. If we now apply the 
Beukers-Kolk-Calabi change of variables [1] 

sin Mi sinw 2 _ sinw 2 n-i _ sinw 2n , Q1fV . 

Xx — , X<i — , • • • , X2n— 1 — j ^2n — (o.lU) 

COS U2 COSW3 ""' COS1i 2n COS Mi 

to the integral ( I3.9p . after some simple manipulations we get 



C(2n+1) 



22n+l _ I \ 2 J 



2n 



In 



7T 

tan ( Mi — 



dux ■ ■ ■ du2n- (3.11) 



A 2 „ 



Using the kernel K2> n (u,v), we can reduce the evaluation of ( 13. lip to the 
evaluation of the following one-dimensional integral: 



C(2n + i; 



2vr 



2n 



2 2n+l _ I 



In 



71 

tan ( —u 
2 



K 2n (u,u) du. 



(3.12) 



But 



In 



tan ( ^(1 - u 



In 



cot ( —u 



— In 



7T 

tan ( — u 
2 



which enables us to rewrite (13 . 1 2[) as 



C(2n+r 



7T 



2d 







x / In 


tan(| U ) 


Jo 





2 2n+l _ ! 

[i^ 2n (u, u) - if 2 „(l - u, 1 - «)] du. (3.13) 
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However, from ( 13. ip and (13. 2p we have (recall that £ , 2n _i(l— u) = —E 2n -i(u)) 

2 2n ~ 1 d 

K 2n (u, U) - K 2n (l It, 1 u) = (-1)™ — -yy — E 2n ( U ), 

and the straightforward integration by parts in (13.131) yields finally the result 

« 2 " +i )- 4 [1 ^t ) ; i (2 „ ) , ^ 

This is exactly the integral representation for ((2n + 1) found in [3]. 



4 Concluding remarks 

Finally, let us comment on the origin of the highly non-trivial Beukers-Kolk- 
Calabi change of variables. Using the hyperbolic version of it 

sinhwi sinhu 2 sinhw n _i sinhii n 

X\ = - , x 2 - , . . . , X n —\ - , x n - , l^.ij 

coshti 2 coshii 3 coshw n coshi^ 

we can get [UJ 

on r r nn 

^ = F^T J ' ' J dUl ' ' ' dUn = 2^T Vo1 ™^' ^ 

u„ 

where U n has a complicated amoeba-like shape with narrowing tentacles go- 
ing to infinity. It was shown by Passare [H] that £(2) is really related to the 
amoeba, a fascinating object in complex geometry (|12j. |10j ) introduced in 
the book [6]. We can conjecture that U n is the 2 n -th part of the amoeba of 
a certain Laurent polynomial and the Newton polytope of this polynomial 
is the union, from which the central point (0, . . . , 0) is thrown away, of the 
(suitably rescaled) polytope A n and its 2™ — 1 mirror images under reflec- 
tions Ui — > —U{. For n = 2, the conjecture is valid and the corresponding 
amoeba uncovers the origin of the two-dimensional version of the Beukers- 
Kolk-Calabi change of variables [TT]. We think it is worthwhile to further 
investigate this possible connection between ((n) and amoebas. 
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A Appendix: The proof of (12.11) and (12.21) by 
induction 



If n = 0, ( 12. 21) gives, because of E (x) = 1, 

^(u,!;) = 0(1 -u-v), 

as desired. Let us calculate K 2 (u,v). From ( 1 1.7ft we have 

i 

K 2 (u,v) = J K\(u, U\) K\(ui, v) du\ = 
o 

1-u 

/" „, , , / 1 - u, if u > u, 

= / 0(1 - v - ux) dux = < . .. , 

J [ 1 — U, 11 V < u. 



We can rewrite this as follows 

K 2 (u, v) = (1 - v)0(v - u) + (1 - u)0(u - u), 

and this is exactly what (11.41) gives for n—1 because E\(x) = x — 1/2. 
Now suppose (12 .ip is true for a particular n and let us calculate 

K 2n+1 (u,v)= / K 1 (u,u 1 )K 2n (u 1 ,v)du 1 . 
Jo 

We will use the following property of the Euler polynomials 

^-E n (x) = nS n _i(x), (A.1) 
ax 

which facilitates the calculation of integrals. 
We have 

-l-u 



[l- u-v) I U E 2n ^(^l±^\ d Ul . (A.2) 



Let us make the change of variables 

Ui + V 



2 X - 
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Then the integral equals to 

.(l-u+O/a i f d 



In V dx 



EoJx) dx 



n 



E. 



2n 



1 — U + V 



- E 2n (v 



Therefore, 



l-u 



E 2 n-1 



U\ + V 



9(1 -u-v) 



n 



E. 



2n 



6{ui — v ) du\ 
1 — U + V 



Eojv) 



(A.3) 



Analogously, using the substitution = x and E 2n (0) = 0, we get 



l-u 



E 2n -1 



U\ — V 



1 — u — U 



(A.4) 



Further we have 



l-u 



Em-\ 



Ui + V 



e(l-u-v) J E 2n ^ ( 



9(v — U\)du\ 

Ui+V 



+9(u+v-l) J E 2n -i 
9(1 - u-v 



l-u 



2 

Ui+V 



du\ 
du\ 



n 



EoJv) - E 2n ( ? 



6(u + v-l) 



2/i 



1 — U -\- V 



— E' 



2n 



(A.5) 



and analogously 

-l-u 



E 2n -i 

(u + v -1) 
n 



v — U\ 



9(v — «i) d«i 



0(1 -u-v) 



/7 



2// 



2?/ 



2// 



It + t> — 1 



(A.6) 



The equations (121)1 . itOj) . fjATijl . (fA~5l) and (TA~6]) indicate that if 2 n+iO,w) 
is given by (12. 2p . as desired. 
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Now let us show that the validity of (12. 2p implies that 
K 2n +2(u,v) = / Ki(u,Ui)K 2n+ i(u 1 ,v)dui 



is given by the formula (12. ip . 
We have 

E 2n 



1 — U\ + V 



6{v-u) / E- 



l-v 



->2n 



6(1 — u\ — v) du\ 
1 — U\ + v " 



+6(u - v) 
2 . 



l-u 



E- 



2/i 



1 — U\ + V 



du\ 



2n + l 
2 

2n + 1 



0(u - w) 
0(u-v) 



-^2n+l(^) — ^2n+l 



1 + V 



Ek, 



2n+l 



U + V 



-Ek, 



2n+l 



1+U 



and 



r-(^F) 

2 



0(1 — Ml — f ) dwi 



2n + 1 
2 



0(v-u) 



-^2n+l (0) — E2n+1 



l-V 



2n + l 



6{u - v) 



E 2 n+l 



U — V 



E211+I 



1 -V 



as well as 



l-u 



E-2n 

2 



1 — Ui + v 



2n + l 



6(v - u) 



E2n+1 



6{u\ + v — 1) dui 
u + v 



E 2 n+l(v) 



and finally 



1— u 



2n 



U\ + v — 1 



0(«i + y — 1) c?Ui 



2n + 1 



0(u-u) 



-E2n+1 

9 



v — u 



-E 2n+1 {0) 



In light of f lA~8]) . (EO]l . fOTTOj) and f lATT]) . the integral (EOjl takes the 

form of (12. ip . with n — > n + 1, because 



which follows from £^(1 — x) = (—l) n E n (x). 
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